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Based on Hijing-1.37 simulations, event-wise mean-pt fluc- 
tuations in Au-Au collisions at -^sjvjv = 200 GeV are studied 
with graphical and numerical methods. This study shows that 
jets/minijets could represent a substantial source of mean-pt 
fluctuations, and jet quenching could reduce these fluctua- 
tions significantly. Mean-pt fluctuation measurements are a 
promising probe of correlation structure produced in the early 
stages of central Au-Au collisions at RHIC energies. 
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I. INTRODUCTION 

QCD lattice gauge calculations [1] suggest that a 
phase transition from hadronic nuclear matter to a color- 
deconfined medium or quark-gluon plasma (QGP) may 
occur in relativistic heavy ion collisions at RHIC. Mea- 
surements of multiplicity and transverse-momentum fluc- 
tuations have been proposed to establish the existence 
[2,3] and study the properties [4-8] of a QGP. It was 
suggested that in case of a second-order phase transition 
fluctuations may be reduced dramatically [4-7] , whereas 
for a flrst-order phase transition large fluctuations may 
occur due to localized QGP droplet formation [8] , or more 
gradual density or temperature fluctuations may occur. 
The past two years have seen signiflcant progress in the 
number and quality of preliminary fluctuation and corre- 
lation measurements at both SPS [9,10] and RHIC [11]. 
Further efforts are now required to draw quantitative con- 
clusions about the phase transition and the physical char- 
acteristics of the QGP. 

In this Letter we consider the possibility that a major 
source of event-wise mean-pt or (pt) fluctuations in rela- 
tivistic heavy ion collisions is jet/minijet production, and 
that these fluctuations can therefore be used to study the 
early collision process and the properties of the QGP, in 
a manner complementary to studies of inclusive event 
properties such as transverse momentum distributions 
[12]. 
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According to [13] minijets (semi-hard parton scatter- 
ing with momentum transfers of a few GeV/c) should be 
copiously produced in the initial stage of ultra-relativistic 
heavy ion collisions at RHIC energies. A QGP formed in 
the early stage of such collisions could act as a dissipa- 
tive medium, causing jets/minijets to lose energy through 
induced gluon radiation [12], a process conventionally 
called jet quenching in the case of higher-pt partons. The 
properties of the dissipative medium should determine 
the specific energy loss of jets and minijets [14]. We 
explore the role of fluctuation measurements in under- 
standing the dissipation properties of a color-deconfined 
medium, using Hijing-1.37 [15] to study jet-quenching 
and jet/minijet contributions to (pt) fluctuations. 



II. HIJING EVENT TYPES 

Three sets of Hijing-1.37 events were generated for this 
study: Hijing default with jet quenching on, with quench- 
ing switched off and with jet/minijet production switched 
off. Because jet quenching effects should be most sig- 
nificant in central collisions this study was applied to 
Hijing-1.37 events with impact parameter less than 3 fm, 
equivalent to the top 5% most central collisions. 

Fluctuation analysis was applied to all generated 
charged particles subject to the conditions that the pseu- 
dorapidity and transverse momentum satisfied \ri\ < 1 
and 0.15 GeV/c < pt <2 GeV/c. The total number 
of events, the mean pt and variance crp^ of the charged- 
particle inclusive pt distributions, and the mean h and 
variance a„ of the charged-particle multiplicity distribu- 
tions for the three types of central Hijing events are listed 
in Table I. We analyze (pt) fluctuations for these three 
event classes using graphical and numerical methods. 



III. GRAPHICAL ANALYSIS OF 
FLUCTUATIONS 

Graphical analysis of (pt) fluctuations directly com- 
pares the data distribution with a reference. The gamma 
distribution models statistically-independent particle 
emission [16]. Any difference between data and reference 
would imply that emitted particles are correlated and/or 
event-wise fluctuations in collision dynamics (such as 
minijets) are present. Effects of fluctuating event mul- 
tiplicity in this method require folding the gamma dis- 
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tribution with a negative binomial distribution [16]. The 
reference distribution is formulated as: 

f{{Pt)) = ^ fNBD{n,l/k,n)fr{{pt),np,nb) (1) 

rtmin 

where fNBDinjl/kjfi) and fT{{pt),np,nb) are the nega- 
tive binomial distribution and n-fold gamma distribution 
respectively 

M{p,},np,ni} = p^(nl.{p,))"'-'exp-"'('"l 

The parameters in Eq. (2) are defined as p = Pt/cr|^, 
b = Pt/'^p^ and n/k = a^/n — 1, where pt and cr|^ are the 
mean and variance respectively of the inclusive pt distri- 
bution. This model is related to the central limit theorem 
result that the cumulants of a distribution of independent 
n-samplc means from a fixed parent distribution are re- 
lated to the cumulants of the parent by inverse powers of 
the sample number. The gamma distribution cumulants 
satisfy this relation, and the inclusive pt distribution is 
approximately described by a gamma distribution with 
n = 2. The relations above then follow. The graphical 
comparison is important as a qualitative indicator of ex- 
cess variance and of significant deviations from central 
limit behavior in cumulants higher than the second. 

A. Basic histograms 

Distributions of (pt) for each of three Hijing event types 
are shown in Fig. 1 as histograms together with reference 
distributions based on Eq. (1) plotted as curves. The 
reference distributions arc normalized to have the same 
peak values as the corresponding data distributions to 
emphasize the width comparison. It can be seen from 
Fig. 1 that the Hijing data histograms deviate substan- 
tially from the reference distributions, indicating that 
nonstatistical fluctuations arc significant. Fig. 1 indicates 
that excess fluctuations in events with jets switched off 
are substantially less than those for the other two event 
classes. 

Comparison of the top and bottom panels of Fig. 1 
indicates qualitatively that jet production plays a signif- 
icant role in {pt) fluctuations in the Hijing model. To 
investigate the effect of jet quenching we can compare 
nonstatistical fluctuations in top and middle panels of 
Fig. 1. However, with these histograms it is not possi- 
ble to detect visually a signiflcant difference in the mag- 
nitude of nonstatistical {pt) fluctuations. We require a 
more differential approach. 
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FIG. 1. Distributions of {pt} for Au-Au collisions at -^/snn 
= 200 GeV from Hijing- 1.37 (histograms) compared with ref- 
erence distributions assuming statistically-independent parti- 
cle emission and calculated using Eq. (1) (curves). The top 
panel shows results for jets-on quench-ofF events, the middle 
panel for jets-on quench-on events and the bottom panel is 
for jets-off events. 

B. Difference histograms 

The graphical analysis can be made more sensitive by 
plotting the difference between data histogram and ref- 
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erence distribution, facilitating quantitative comparisons 
between Monte Carlo configurations and references. The 
normalized event-number difference is defined as 

CAT-/ /aT ^data ^reference /o\ 

5NI^/N = (3) 

V^^data 

The independent random variable \/fi{{pt) — 'Pt)/<^pt 
normalizes the distribution width relative to the inclusive 
width according to a central limit expectation. The result 
is a universal plotting format; the horizontal and vertical 
axis units are standard deviations. 
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FIG. 2. Normalized event-number difference SN/^/N vs 
V^{{Pt) — Pt)/o'pt for three types of Hijing-1.37 events, for 
Au-Au collisions at ^/snn = 200 GeV. Solid line shows results 
for jets-on quench-off events, dashed line for jets-on quench-on 
events and dotted line for jets-off events. 



tions. The linear variance comparison measure Acp^ [17] 
was developed to isolate pure (pt) fluctuations from ef- 
fects of multiplicity fluctuations which may appear as a 
systematic bias in some statistical measures. 

Several quantities have been proposed to measure 
event- wise (pt) fluctuations in relativistic heavy ion colli- 
sions: fp, [6], Fp, [11], [18], and 

^pt .dynamical 

[19]. A 

detailed study of is given in [17]. In the present anal- 
ysis the closely- related Adp^ is applied to event- wise (pt) 
fluctuations in Hijing-1.37 events. This difference factor 
for (pt) fluctuations is defined as 

Aap, EE {niipt) - ptf - crl}/2ap, (4) 

where the overline means average over all events in a data 
set. Values of AcTp^ for the three Hijing event types are 
included in Table I. Consistent with Fig. 1 the results 
show very significant contributions from jets, and they 
also show that fluctuations are significantly reduced when 
jet quenching is applied, consistent with Fig. 2. 

As expected, numerical analysis of (pt) fiuctuations 
reveals more precisely than graphical analysis that 
jet/minijet production in the Hijing model contributes 
significantly to nonstatistical event-wise (pt) fluctuations. 
The increase in fluctuations for jets on compared to jets 
off is a factor 3-5. The effect of jet quenching is to re- 
duced AcTpj by about 30% relative to no quenching, 20 
times the statistical error. These results suggest that 
properties of the QGP or dissipative medium formed in 
the initial stage of heavy ion collisions may be studied 
precisely with numerical analysis of event-wise (pt) fluc- 
tuations. 



Figure 2 shows results for the three Hijing-1.37 configu- 
rations. Significant deviations from zero reveal nonstatis- 
tical fluctuations in all three event types. The dramatic 
increase from jet-off events (lowest curve) to events with 
jets on is consistent with the qualitative difference in Fig. 
1. A significant difference with and without jet quenching 
(upper two curves) is also now apparent. This difference 
indicates that according to the Hijing model transverse 
momentum fluctuations are sensitive to jet quenching, 
which may in turn reflect the properties of a dense col- 
ored medium or QGP formed in relativistic heavy ion 
collisions. We conclude that while the original graphical 
method is an indicator for fluctuations associated with jet 
production, the difference histograms are more sensitive 
to the effects of jet quenching. 

IV. NUMERICAL ANALYSIS OF 
FLUCTUATIONS 

Numerical analysis of nonstatistical fluctuations re- 
quires comparison of a variance from data with the vari- 
ance of a reference representing only statistical fluctua- 



V. CENTRALITY DEPENDENCE 

However, we observe an excess of binary-collision scal- 
ing in Hijing compared to RHIC data which affects the 
interpretation of this fluctuation analysis. Hijing particle 
production increases with number of binary collisions 

TABLE L The total number of events, mean pt and vari- 
ance (Tpj of inclusive charged-particle transverse-momentum 
distributions, the mean n and variance a„ of charged-particle 
multiplicity distributions, and the event-wise mean-pt excess 
fluctuations Aap^ for Hijing-1.37 events — 5% central (b < 3 
fm) Au-Au collisions at y^ijvjv = 200 GeV, |r;| < 1. Errors 
are statistical only. 



event class 


no jets 


jets - quench on 


jets - quench off 


event total 


134868 


121671 


146241 


n 


431 


1387 


1159 




30 


130 


100 


Pt(MeV/c) 


424 


478 


496 


<7p,(MeV/c) 


227 


283 


307 


Acrp,(MeV/c) 


7.2 ±0.3 


24.9 ±0.5 


34.5 ±0.5 
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(ccntrality) at twice the rate observed in RHIC Au- 
Au collisions, as indicated by the multiplicity entries in 
Table I. The jets-off multiplicity reflects soft processes, 
whereas for the same number of participant pairs the jets- 
on multiplicities, which include hard-component scaling 
with binary collisions, represent a 2.5-3-fold increase over 
participant scaling for central events. For RHIC data we 
observe a corresponding 1.5-fold increase [20]. The frac- 
tion of particle production from hard processes in Hijing 
events is thus about twice that observed for RHIC events. 
We find that {pt) fluctuations in Hijing (jets on) are dom- 
inated by the hard component, by a factor 5. It may be 
advisable therefore to reduce by a factor of two the (pt) 
fluctuations attributed by this study to jets/minijets in 
Hijing in comparisons with RHIC collisions. 

VI. HADRONIC RESCATTERING 

Hadron rcscattering between chemical and kinetic de- 
coupling might attenuate nonstatistical fluctuations gen- 
erated in the early stage of the collision, which could then 
be confused with effects of dissipation in the prehadronic 
medium. Recent Monte Carlo calculations suggest that 
hadronic rescattering may be negligible [3,21]: the time 
for rescattering may be so brief [22,23] that fluctuation 
attenuation is not signiflcant. Nonstatistical {pt) fluctu- 
ations due to jets/minijets would thus survive hadronic 
rescattering unscathed. This is an important point for 
interpreting attenuation of fluctuations observed in the 
final state. Reduced {pt) fluctuations as reported in this 
Letter would then be an unambiguous indication of pre- 
hadronic dissipation (minijet/jet quenching). 

VII. CONCLUSIONS 

In summary, event- wise {pi) fluctuations in central Au- 
Au collisions at ^Jsnn = 200 GeV have been studied 
with Hijing-1.37 using graphical and numerical methods. 
The study shows that jet/minijct production in the ini- 
tial stage of heavy ion collisions may contribute substan- 
tially to nonstatistical {pt) fluctuations. We have also 
shown that in the Hijing model (pt) fluctuations are sen- 
sitive to jet quenching, thus suggesting that fluctuation 
analysis could serve as a probe of dissipation properties 
of a dense colored medium formed in the initial stages 
of ultrarelativistic heavy ion collisions. Studies of fluc- 
tuation centrality dependence in models and data may 
help clarify questions on hadronic dissipation and Hijing 
binary-collisions scaling raised by this study. 

We appreciate helpful discussions with J.G. Reid (Uni- 
versity of Washington) . This work was supported in part 
by USDOE contract DE-FG03-97ER41020. 
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